Lattice distortions, incommensurability, and stripes in the electronic model for 

high-T c cuprates 
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Striped superconductivity (SC) with lattice distortions is investigated based on the three-band 
Hubbard model for high-T c cuprates. A stable inhomogeneous striped state is determined in the 
low-temperature tetragonal (LTT) phase with lattice distortions using a quantum variational Monte 
Carlo method. The ground state has vertical or horizontal hole-rich arrays coexisting with incom- 
mensurate magnetism and SC induced by several percents of lattice deformations. The SC order 
parameter oscillates according to the inhomogeneity in the antiferromagnetic background with its 
maximums in the hole-rich regions, and the SC condensation energy is reduced as the doping rate 
decreases. 

PACS numbers: 78.20.-e, 78.30.-j, 74.76.Bz 
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Over the last decade the oxide high-T^superconduc- 
tors have been investigated intensively, [lj The mech- 
anism of superconductivity (SC) has been extensively 
studied using various two-dimensional (2D) models of 
electronic interactions. The 2D three-band Hubbard 
model is the simple and most fundamental model among 
such models. The 2D one-band Hubbard model is re- 
garded as the simplified model of the three-band model. 
Studies of these models over the last decade indicated 
that the d-wave SC is induced from the electronic repul- 
sive interactional^, 0,0 S0 ; significantly it has been 



0.004 



LU 

< 



0.003 



0.002 



0.001 




0.005 



0.01 



1/N. 



FIG. 1: SC condensation energy per site as a function of 
1/N a in t units where t w t p d/3 and -/V a is the number of 
atoms. Squares are for 8 w 0.2, t pp /t p d = 0.4 and Ud/t v d = 8 
for the three-band model on square lattices. Circles are at S — 
1/8 coexisting with stripes for t pp /t p d = 0.4 and Ud/t p d = 8 
on rectangular lattices 32 x 8, 24 x 6, 16 x 8 and 16 x 4. 
Triangles are for the single-band Hubbard model; S — 0.86 
and t' — —0.2 and [7 = 8 for solid symbols and 5 = 0.84 
and t' = —0.15 for open symbols (energy unit is £).[jj The 
diamond shows the value indicated from experiments. 



shown that the SC condensation energy and the mag- 
nitude of order parameter are in reasonable agreement 
with the experimental results in the optimally doped 
case, jsl tlOl tllj 

The SC condensation energy obtained by the vari- 
ational Monte Carlo method (VMC) is estimated as 
E con d — 0.00117t = 0.59meV per site in the opti- 
mally doped case for the single-band Hubbard model 
in the bulk limit. 0, We must note that E con d is 
given as 0.17 ~ 0.26meV by specific heat data[l2l 
lia ] and 0.26meV by the critical magnetic field value 
-ffc/ 87r Q- The agreement of the VMC value with 
the experimental estimation is quite significant and sup- 
ports the calculations. The VMC method can be re- 
garded as an approximation to Quantum Monte Carlo 
calculations. |l5l fill The SC order parameter A s deter- 
mined from a minimum of the energy is of the order of 
0.01-0.015=15meV-20meV at hole density S ~ 0.2. \H\ 

The interplay between magnetism and superconduc- 
tivity is suggested in the underdoped region. The re- 
duction of T c in this region remains unresolved and may 
be related to magnetism. An existence of incommensu- 
rate correlations with modulation vectors given by Q s = 
(tt ± 27n5,7r) and Q c = (±4tt£,0) (or Q s = (7r,7r ± 2tt5) 
and Q c = (0, ±47r<5)) has been sug gested for the hole- 
doping rate M, [H HI M, ll E3 The linear 
doping dependence of incommensurability in the under- 
doped region supports a striped structure and suggests 
a relationship between magnetism and SC.0| A rela- 
tionship between the SDW, CDW orders and a crys- 
tal structure is also suggested in intensive studies by 
the neutron-scattering measurements; ^3, 0> HE US m 
particular in the low-temperature tetragonal (LTT) and 
low-temperature less-orthorhombic (LTLO) phases, the 
CDW order is stabilized, [27j while no well defined in- 
commensurate CDW peaks were observed for the or- 
thorhombic systems. [H H|| In the elastic and inelastic 
neutron scattering experiments with I^-^Sr^CuO^ the 
incommensurate magnetic scattering spots around (7r,7r) 
have been observed in the SC phase in the range of 
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In this paper an incommensurate striped SC under lat- 
tice distortions is investigated based on the three-band 
model. In the description of SC in the underdoped re- 
gion it is of prime importance to investigate the effect 
of inhomogeneity. In this paper we take into account all 
of the inhomogeneity, lattice instability and anisotropic 
pairing to clarify the ground state in the underdoped re- 
gion of high-T c cuprates. The Hamiltonian is given by 

H = H pd + V + H el where 

icr icr 

icr 

- Pi-x/2,a - Pi-y/2,<r) + h.C.] 

+ tpp ^(1 + V. t )[pl + ~ /2 a p i+s ./ 2 ,a ~ p\ + y/2,,jPi-x/2,a 
icr 

~ pl-y/2,aPi+x/2,a + Pi_ v / 2j(7 Pi-x/2,a + h.C.]. 
+ ipd^^[Uixd\ a Pi+x /2,a ~ u i,-xd\ a Pi- x /2,a 



+ u iydl a Pi+y/2,a - Ui-yd\ a Pi-yl2,a + h.C. 



y-direction, respectively, p\±%i 2 a an( ^Pi±x/2,a denote the 
operators for the p electrons at the site Ri ± x/2, and in 
a similar way p\ ±y / 2 a an< ^ Pi±y/2,a are defined. Ud de- 
notes the strength of Coulomb interaction between the 



d electrons. 



t! ifi 



and vt represent the variations of the 



transfer energy t p d and t pp , respectively. The number of 
cells which consist of d, p x and p y orbitals is denoted as 
N. 

The wave function with the inhomogeneous spin 
structure is made from solutions of the Hartree-Fock 
Hamiltonian given as H tria i = H^ d + J2iA^ n di ~ 
(7{—l) Xi+Vi mi]d\ cr di ai where we have variational param- 
eters €p and id in H® d . In this paper Sridi and 
are assumed to have the form 



J2j a/cosh((x l 



")), and m l = A inc J7 - tanh((a 



Xj tr )), with parameters a and Aj nc where xf r denote 
the position of a stripe. The inclusion of stripe order pa- 
rameters considerably improves the ground-state energy. 
In small clusters the deviation of the energy of striped 
state from the exact value is within several percents for 
the Hubbard model. |33| 

The wave function is constructed from the solution of 
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where H e i denotes the lattice elastic energy given 



by H e 



(K pp /2)J2i^ v i where K p d and K pp denote the elastic 
constants, x and y represent unit vectors in the x- and 
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FIG. 2: SC condensation energy per site vs the hole density 
in tpd units, where the parameters are t pp = 0.4 and Ud ~ 8. 
Solid circles and open circles indicate the SC condensation 
energy for the uniform SC and striped SC, respectively. The 
lines are fitted by parabola. Squares are obtained for the 
single-band Hubbard model with the next-nearest transfer 
t' = -0.2 on 12 x 12 lattice. ^ 



FIG. 3: Lattice structures in the LTT phase (a) and LTO 
phase (b). The symbol "+" means that the oxygen atoms 
move upward and instead "-" oxygen atoms move downward. 
"O" denote the oxygen atom. 
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Bogoliubov-de Gennes equation given by 

J2(H m u$ + F ijV j) = E*u}, (3) 

3 

J2( F ^4 - H MV$) = E x vl (4) 
i 

where (Hij a ) and (Fij) are 3iV x 37V matrices including 
the terms for d, p x and p y orbitals. The Bogoliubov 
operators are written in the form 

ax=J2( u i a H+ v i a li) (E X >0), (5) 

i 

«a = E( u ^t+ w H) (^<o). (6) 

i 

ai a denotes di a , Pi +X /2<r or Pi+#/2o- corres P on ding to the 
components of u x and . 

Then the wave function is written as [34ll35l l3^| 

l/> = P G PN c Y[a x a[\0) 

A 

cx PGiY^iU-'V)^]^/^). (7) 

ij 

U and V are matrices defined by (V)\j = v x and (U)\j — 
Uj. Pq is the Gutzwiller operator. The spin modula- 
tion potential is contained in (Hij a ) and the SC order 
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FIG. 4: Energy gain AE = E(u = 0) - E(u) (u = 5t pd /t pd ) 
per site as a function of transfer deformation u in t pd units. 
The parameters are t pp = 0.4 and Ud = 8 for 16 x 4 lattice. 
The hole-rich stripes are in the {/-direction. The energy gains 
for (A) Uiy = (triangles), (B) ua = 0, (open circles) and (C) 
(solid circles) are shown. The elastic energy Ku 2 /2 is shown 
by the dashed line (for K — 5 and K — 10) . The summations 
of AE — E(u = 0) — E(u) and the elastic energy per site 
are also shown for the case (C). The Monte Carlo statistical 
errors are smaller than the size of symbols. 



parameters Ay are included in (Fij). We assume the 
following spatial variation for the SC order parameters 
in the (i-electron part: A»,i+s = A s cos(Q x (xi + x/2)), 
&i,i+y — —A s cos(Q x Xi), where Q x = 2t:S (S is the hole 
density). The SC order parameter oscillates according to 
the spin and charge distributions so that the amplitude 
has a maximum in the hole-rich region and is suppressed 
in the hole-poor region. The energy-expectation value 
is calculated using the Monte Carlo algorithm jrl ITU Hi} : 

Here we show the results in the case without lattice 
distortions. It has been shown that the striped state is 
more stable than the uniform SDW state for small hole 
doping. 32] In Fig.l the size dependence of the SC con- 
densation energy is shown for the uniform SC in the over- 
doped region and the striped SC in the underdoped re- 
gion with the results obtained for the one-band Hubbard 
model for comparison. The parameters are t pp = 0.4 and 
Ud = 8 in t p d units. The squares in Fig.l indicate the 
SC condensation energy of pure ei-wave state at S w 0.2, 
while the circles are for SC coexisting with stripes at 
S = 1/8 for Q x — 7r/4 evaluated on rectangular lattices 
32 x 8, 24 x 6, 16 x 8 and 16 x 4. In both cases the energy 
obtained through an extrapolation is of the same order 
as experimental values. 

E cond P3 0.000Ut p d w 0.2meV, (8) 

where we have assigned t p d « 1.5eV.[32j The data in 
Fig. 2 show the SC condensation energy as a function of 
the hole density. The SC condensation energy per site for 
the striped SC is reduced as the hole density decreases, 
while that for pure d-wave SC remains finite even near 
half-filling. This suggests that an origin of the decrease 
of T c in the underdoped region lies in the reduction of 
hole-rich domain where the SC order parameter has finite 
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FIG. 5: Schematic structure of lattice distortions and stripes 
where the hole-rich arrays are perpendicular to the tilting 
axis. We call this state the mixed LTT-HTT phase. The 
shaded square represents distorted CuO unit cell. 
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amplitude. 

Now let us consider the effect of lattice distortion on 
stripes. In the LTT phase stabilized at low tempera- 
tures near 1/8-hole filling, [3^] the distortions of the CuO 
square occur in the manner shown in Fig. 3. The LTT 
phase has a 'tilting axis' on which the copper and oxy- 
gen atoms never move even in the distorted state. The 
vertical or horizontal stripes can coexist with the lattice 
distortions in the LTT phase, being parallel or perpen- 
dicular to the tilting axis. 

The structural transition from low-temperature or- 
thorhombic (LTO) to LTT phases occurs in LaBaSrCuO 
and LaNdSrCuO systems. |24j It is not clear a priori what 
structure is stabilized due to the lattice deformation. We 
consider the following cases assuming that the stripes 
are in the y-direction: 

(A) u ix = u, Uiy = 0, Vi = 0, 

(B) ua = 0, u ly = u 1 v l = 0, 

(C) Ui$ = u,Uiy = ucos(2Q x x i ),u l = ucos(2Q x Xi), 

where Q x = 2ttS and u is the amplitude of defor- 
mation: — Ui t —p,. u = corresponds to the LTO 
sturcture, and the anisotropy in t p d indicates a transi- 
tion to the LTT phase. t p d increases along the tilt axis 
compared to the LTO phase. The case (C) corresponds 
to the structure of mixed LTT-HTT phase. The energy 
gain per site defined as AE/N = (E(u = 0) - E(u))/N 
is presented in Fig. 4 as a function of u in t p d units. The 
energy in the case (B) is lower than that in the case 
(A) indicating that the stripes are parallel to the tilting 
axis under the rigid LTT structure. The energy in the 
cse (B) is lower than that in the case (A) indicating 
that the stripes are parallel to the tilting axis under 
the rigid LTT structure. We simply assume the same 
elastic energy cost for these types of rotations. The cost 



of energy due to lattice distortions is assumed to be 
given by (K/2)u 2 for the constant K. K is estimated as 
follows. According to Harrison's rule,^ tpd is expected 
to vary as d~ n with n « 7/2, d being the Cu-0 bond 
length. Since St p d/t p d = —nSd/d, the elastic energy is 
estimated as 

E el = \c{2df{^ f = \c{2df^ ee | u 2 . (9) 

The constant C is estimated as C rj 1.7 x 10 12 dyn/cra 2 = 
1.7 eV/1 3 . 41] Since d ss 2 A, K is of the order of lOeV: 
K « 8.9eV. We point out that E e i has possibly a linear 
term in u. The 9 2 term, if present in E e i, is proportional 
to u since u ~ 1 — 003(6*) ~ 8 2 /2 where 8 is the tilt 
angle. The presence of linear term may lead to a first 
order transition. As shown in Fig. 4 the striped state is 
more stabilized in the LTT phase. We show schemati- 
cally the stable striped state in the LTT phase in Fig. 5 
obtained from our VMC evaluations, where the shaded 
square represents the tilted CuO unit cell rotating around 
the tilting axis. The LTT-HTT state is more stabilized 
due to the kinetic energy gain coming from the softening 
of tilt angles. 

In this paper we have investigated the inhomogeneous 
ground state with the lattice distortions based on the 
three-band model of high-T c cuprates using the varia- 
tional Monte Carlo method. The SC condensation energy 
decreases as the doping ratio decreases, which is due to 
the reduction of SC domain in the hole-rich region. The 
stable striped state has hole-rich arrays being perpendic- 
ular to the tilting axis of the lattice distortions in the 
LTT phase as shown in Fig. 5, which can be regarded as 
the LTT-HTT mixed phase. We thank H. Oyanagi for 
valuable discussions. 
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